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We review the level spectroscopy, which is a powerful method of analyzing the numerical 
data with respect to the Berezinskii-Kosterlitz-Thouless quantum phase transition in one 
dimension. We focus on its physical meaning and also its application to the magnetization 
plateau problems. 

§1. Introduction 

The Berezinskii-Kosterlitz-Thouless ^ " 4 ) (BKT, or simply KT) quantum phase 
transition in one-dimensional systems is the transition between the gapless state 
and the gapful state. As is well known, the critical behavior of the BKT transition 
is highly singular and sometimes called "pathological". In finite systems, this high 
singularity appears as severe slowly-converging logarithmic size corrections in various 
physical quantities. Thus it is very difficult to determine the BKT quantum phase 
transition point from the numerical data if we use the conventional methods. 

We note that there are two types in the BKT quantum phase transition. One 
is the doubly degenerate type, of which examples are (a) the transition between the 
spin-fluid (SF) and Neel states in the 5 = 1/2 XXZ spin chain, and (b) the SF- 
dimer transition in the 5 = 1/2 Heisenberg spin chain with next-nearest-neighbor 
interactions. 5 ) " 8 ) In these cases, the gapped states (the Neel state of (a) and the 
dimer state of (b) are doubly degenerate. The mechanism of the transition of this 
type is the spontaneous symmetry breaking. Another type is the non-degenerate 
type. The examples are (c) the SF-Haldane transition in 5 = 1 XXZ spin chain, 
and (d) the SF-dimer transition in 5 = 1/2 bond-alternating XXZ spin chain. 9 )' 10 ) 
In the non-degenerate case, the gapped states (the Haldane state in (c) and the dimer 
state in (d)) are unique and non-degenerate. The nature of this type of transition 
is that the gap- formation mechanism (for instance, the bond alternation in (d)) in 
the Hamiltonian is renormalized to zero in the sense of the renormalization group 
treatment due to the strong quantum fluctuations. 

In this paper, we review the level spectroscopy 11 '' 12 ) (LS) by use of which we 
can determine the BKT critical point very accurately (typical accuracy is 1CP 3 or 
better) from the numerical diagonalization data, overcoming the above-mentioned 
difficulties. We mainly focus on the physical meaning 13 ) of the LS for the doubly 
degenerate type in the one-dimensional quantum spin models. Also we explain how 
to apply the level spectroscopy to the magnetization plateau problems. 
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§2. Physical Meaning of the Level Spectroscopy 

To distinguish the gapless spin-fluid state and the gapful state, the most funda- 
mental quantity is the excitation gap for infinite systems. If we know the analytically 
exact solution, we can easily distinguish the gapless and gapful states. In usual cases, 
however, we do not know the exact solution, so that we have to rely on the numerical 
methods, for instance the numerical diagonalization. Although we may extrapolate 
the gap data to N — > oo (N is the system size), the extrapolated data contains nu- 
merical errors, which brings about the difficulty in judging whether the extrapolated 
value is zero or finite. Further, the extrapolation from the N < £ (£ is the correlation 
length) data is unreliable. If we use this type of extrapolation, the apparent BKT 
critical point will be the point of £ ~ iVo, where Nq is the typical system size of the 
diagonalization data. Thus this method results in the wider gapful region than the 
true one. 

The shortcoming of the extrapolation method lies in the zero-or-finite judgment. 
If we find two quantities which cross with each other at the BKT critical point 
as functions of quantum parameters in the Hamiltonian, we can easily know the 
candidate for the critical point. This is the fundamental idea of the LS. 
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Fig. 1. The finite effects are schematicaiiy shown for the (a) doubly degenerate ground state (gap- 
less) and (b) the spin-fluid ground state (gapless). GS is the abbreviation of "ground state". 
The shadowed region represents the continuous spectrum. 



In usual cases, the ground state of the antiferromagnetic spin model is unique 
(non-degenerate) for finite systems, even when the degenerate ground state is realized 
for infinite systems. The exceptions are, for example, Ising model and the Majumdar- 
Ghosh model. Then, how does the doubly degenerate ground state take place? The 
mechanism is as follows. A low-lying excitation of the finite system asymptotically 
degenerate to the ground state as N — > oo. In this limit, a recombination of these 
two states occurs, which results in the realization of the doubly degenerate ground 
states with spontaneous symmetry breaking. In other words, the double degeneracy 
in infinite systems is lifted by the perturbation of 'finiteness' as is schematically 
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shown in Fig. 1(a). Then the lowest excitation in the gapful region should be the 
broken half of the doubly degenerate states denoted by A in Fig. 1(a). The property 
of the A state depends on what kind of doubly degenerate state is realized in infinite 
systems. 

On the other hand, the gapless spectrum in the spin fluid case becomes discrete 
for finite systems. The lowest excitation should be the spin-wave excitation denoted 
by B in Fig. 1(b). 

Then the properties of the lowest excitations for finite systems in case of the SF 
and gapful cases are quite different from each other. Therefore the BKT critical point 
can be obtained from the crossing of these two excitations as functions of quantum 
parameters in the Hamiltonian. Thus our method is named level crossing method, 
or more sophisticatedly level spectroscopy. 



§3. An Example: 5 = 1/2 XXZ Spin Chain with 
Next-Nearest-Neighbor Interactions 



As an example, let us take up the S =1/2 XXZ spin chain with next-nearest- 
neighbor interactions 5 ) " 8 ) described by 

H = £ {(Sj ■ S j+1 ) A + a(Sj ■ S j+2 ) A } (3-1) 



where 



(S l ■ S m ) A = Sf S* + SfSl + AS?S* 



(3-2) 
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Fig. 2. Phase diagram of the model (3.1). 



For simplicity we restrict ourselves to 
the < A < oo and < a < 1/2 
case. In this parameter region, the 
ground state phase diagram consists of 
three phases, the SF phase, the dimer 
phase and the Neel phase, as is schemat- 
ically shown in Fig. 2. The SF state 
is unique and gapless, whereas the lat- 
ter two states are doubly degenerate and 
gapful. 

We consider the A = 1 case where 
the total spin 5 t ot is a good quantum number, for example. When N = 4n case 
(n is an integer), the ground state for the finite systems is of Stot = 0. Since the 
dimer ground state is also of Stot = 0, the lowest excited state in the dimer region, 
which is nothing but the broken half of the doubly degenerate ground states for 
infinite systems, should be also of Stot = due to the addition rule of the angular 
momentum. On the other hand, the lowest excitation in the SF region should be of 
the spin- wave type with S to t = 1 (triplet due to the SU(2) symmetry). Thus the 
crossing of the S to t = and S to t = 1 excitations is the candidate for the SF-dimer 
critical point, as schematically shown in Fig. 3. 
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Since the crossing point is weakly 
dependent on N, the final result is 
known by taking liniAr^oo a c (N). By 
this procedure, we can obtain a c = 
0.2411 for A = 1. We note that the 
only 4 spin result is a c (N = 4) = 0.25, 
which suggests the effectiveness of the 
LS method. 

When A 1, 5 to t is not a good 
quantum number. In this case, we 
should classify the excitations by dis- 
crete eigenvalues of the symmetry oper- 
ators. They are 5 t z ot , the wave number 




Fig. 3. Schematic behavior of the dimer exci- 
tation with Stot = (solid line) and the 
spin-wave excitation with Stot = 1 (dot- 
ted line). The ground state should be the 
dimer state or the SF state depending on 
whether the former is smaller or larger then 
the latter. 



k, the eigenvalues of the parity operator 
P (Sj <-> Sjv-j+1) and the time reversal 
operator (| f) <-> | J,)). Similar consider- 
ation to the A = 1 case leads to the classification of the excitations as shown in Table 
1. The T operator is valid for the 5 t z ot = case. We can also find the dimer- Neel 
transition point from the crossing of the dimer and Neel excitations. This dimer- Neel 
transition is of the second order with the Gaussian universality in which the critical 
exponents varies continuously along the transition line. 

§4. Application to the Magnetization Plateau Problems 



Table I. Symmetry classification of the ground 
state and important excitations. 
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Since the magnetization plateau is 
related to the field-induced excitation 
gap, we can apply the LS method to 
this problem. There are also two types 
of mechanisms in the magnetization 
plateau, as in case of the zero magneti- 
zation case. The examples of the doubly 
degenerate type are (a) Ms/2 plateau 
in the 5=1/2 Heisenberg chain with 

bond-alternation and the next-nearest-neig hbor interactions, 14) - 15) and (b) (2/3)M s 
plateau in the 5=1/2 distorted diamond chain. 16 )' 17 ) For the non-degenerate case, 
the examples are (c) Mg/3 plateau in the 5 = 1/2 ferromagnetic- ferromagnetic- 
antiferromagnetic spin chain, 18 )" 20 ) and (d) Ms/ plateau in the 5 = 3/2 XX Z spin 
chain with on-site anisotropy. 21 )' 22 ) We note that the model (c) is the first model in 
which the magnetization plateau in quantum spin chain is discussed. 

Because we should compare the lowest energies Eq with the magnetization M = 
Mo, Md ± 1 in the plateau problems at the magnetization Mo, a slight modification 
from the no magnetic field case is necessary. The difference in M corresponds to 
the difference in the number of fermions in the fermion picture. Thus we have to 
take the effect of the chemical potential fi into consideration. Namely, we should 
compare E (M ) and Eq(M + 1) — /x, and also E (M ) and £t)(M — 1) + fx. Since 
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the chemical potential around Mq is given by 

2ju = E (M + 1) - E (M - 1) (4-1) 

we see 

{E (M + 1) - n} - E (M ) = {E (M - 1) + M } - £ (M ) 

= (1/2){£ (M + 1) + £ (M - 1)} - £ (M ) (4-2) 

which plays the role of the spin-wave excitation in the zero magnetic field case. 
The same expression can be obtained by use of the Legendre transformation E — ► 
E-MH. 

Using this method, we have obtained the plateau phase diagrams of the (2/3)Mg 
plateau of the 5 = 1/2 distorted diamond spin chain, 16 )' 17 ) the Ms/2 plateau of the 
5 = 1/2 frustrated two-leg spin ladder, 23 ) the Ms/2 plateau of the 5 = 1/2 two- 
leg spin ladder with 4-spin cyclic exchange interactions, 24 ) the Ms/4 plateau of the 
5 = 1 frustrated two-leg spin ladder. 25 ) " 27 ) Yamamoto, Asano and Ishii 28 ) applied 
this method to the Kondo necklace model with next-nearest-neighbor interactions in 
a successful way. 

§5. Summary 

We have explained the level spectroscopy (LS) especially focusing on its physical 
meaning in the doubly degenerate cases. Not only we can obtain the BKT critical 
point from the level crossing, but also we can check the universality class by the 
careful examination of the combination of the low-lying excitations to eliminate the 
troublesome logarithmic size corrections in the lowest order. 8 )' 29 ) Thus the name 
level spectroscopy may be more appropriate than the name level crossing method for 
our method. 

For the non-degenerate case, of which we have mentioned only briefly, the LS 
method is also established, 30 ) " 34 ) although it is more complicated and the physical 
meaning is not so clear as the degenerate case. The use of the twisted boundary 
condition 30 ^' 31 ) has been proved to be powerful in many cases of the zero-magnetic 
field BKT transitions as well as the plateauful-plateauless transitions. 20 )" 22 )' 35 )" 38 ) 
This twist boundary condition method is also applicable to the transition between 
two plateau states with different plateau formation mechanisms. 21 )' 39 )' 40 ) 

The LS method has been applied not only to the quantum spin models but also 
to the lattice electron models, for instance, t — J and related models, 41 )" 43 ) and the 
extended Hubbard and related models. 43 )" 49 ) 
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